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VIBRATION  OF  COMPOSITE  STRUCTURES 
Charles  W.  Bert 

School  of  Aerospace,  Mechanical  and  Nuclear  Engineering 
The  University  of  Oklahoma 


1. 


INTRODUCTION 


The  purpose  of  this  paper  is  to  review  some  of  the  recent  contributions  to 
the  field  of  vibration  of  composite  structures.  The  topics  of  fatigue  and  frac¬ 
ture,  while  quite  important  in  the  overall  dynamic  performance  of  composite 
structures,  are  not  considered  in  this  survey. 


The  survey  treats  these  topics:  dynamic  stiffness  and  damping  of  composite 
materials,  beams  and  curved  bars  and  rings,  flat  panels,  and  cylindricallv  curved 
panels  and  shells.  It  concludes  with  some  suggestions  for  future  research.. 


2.  DYNAMIC  STIFFNESS  AND  DAMPING. QF  COMPOSITE  MATERIALS 
2 . 1  Stiffness  of  Unidirectional  Fiber-Reinforced  Material 


It  is  well  established  that  fiber-reinforced  composite  materials,  by  virtue 
of  the  geometry  of  the  fibers  acting  as  stiffeners,  behave  macroscopically  like 
homogeneous  but  anisotropic  materials.  The  class  of  elastic  symmetry  depends 
upon  the  geometric  arrangement  of  the  fibers  ( Table  1). 

/f" 

Table  1.  Effect  of  Fiber  Arrangement  on  Elastic 
Symmetry  of  Unidirectional  Composites 


Arrangement  of  Fibers 
in  Cross  Section 


Practical 

Examples  Class  of  Elastic  Symmetry 


Very-small-diameter  fibers 
packed  statistically  iso¬ 
tropically 


Glass/resin,  Transversely  isotropic  (isotropic 

graphite/resin  in  fiber  cross-sectional  plane) 


Hexagonaiiy  packed  fibers  Same 


Same 


Fibers  arranged  in  a  single  3oron/epoxy 
row  or  in  a  rectangular  tape 

array 


Orthotropic  with  planes  of  sym¬ 
metry  intersecting  in  the  fiber 
direction  and  directions  of  rows 
and  columns  of  fibers 


Assuming  linear  stress-strain  behavior,  the  generalized  Hooke's  law  may  be 
written  in  compliance  matrix  form,  using  contracted  notation,  as  follows: 

Y  •  Y:'Y  <u 

where  in  the  three-dimensional  case,  the  strain  matrix  and  stress  macrix  are 
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Here  =  normal  strain,  y^j  =  engineering  shear  strain,  =  normal  stress, 
and  t.-j  =  shear  stress. 


For  the  orthotropic  case,  the  compliance  matrix  takes  the  fora 
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(3) 


For  the  case  of  transverse  isotropy,  with  the  23  plane  as  the  plane  of  iso¬ 
tropy,  array  (3)  takes  the  following  special  fora:  S 1 3  =  S>2»  S33  ■  s22> 

=  2(S^;-S'_2)  and  S 5 5  =  S 5 5  . 


Equation  (1)  can  be  inverted  to  obtain  the  3-0  stiffness  fora  as  follows: 


(4) 


Here,  the  C^_j  are  the  three-dimensional  Cauchy  stiffness  coefficients.  Although 
equation  (1)  is  applicable  to  one-,  two-,  or  three-dimensional  stress  systems, 
equation  (4)  is  applicable  only  in  the  case  of  3-D  stress  systems.  In  the  case 
of  a  relatively  chin  member,  it  is  often  acceptable  to  neglect  the  effect  of 
thickness  normal  stress  (03),  then  equation  (1)  is  inverted  to  obtain  the 
following  expression 


ia . 


(5) 


where  the  Q^j  are  the  plane-stress  reduced  stiffnesses  which  are  related  to  the 
three-dimensional  Cauchy  stiffnesses  as  follows: 

Q. .  =  C. .  -  (C.,C.,/C_.) 
ij  ij  13  jj  :: 

All  of  the  above  elastic  coefficients  (S^j,  C<^, 
the  so-called  engineering  coefficients,  Young's  moduli 
and  shear  moduli  G^j ;  see,  for  example,  [i,2,3j. 

It  is  noted  that  both  the  compliance  [S^]  and  stiffness  ([Cjjj  or  [Q^j]) 
matrices  are  symmetric  as  a  consequence  of  conservation  of  strain  energy,  as  was 
first  shown  very  elegantly  by  George  Green  in  the  early  1800's.  This  reduces  the 
number  of  independent  elastic  coefficients  in  every  case.  It  is  emphasized  that 
in  writing  all  of  the  preceding  equations  the  stresses,  strains,  and  elastic  co¬ 
efficients  are  all  related  to  material-symmetry  axes,  i.e.,  axes  formed  by  the 
intersection  of  mutually  perpendicular  planes  of  symmetry .  In  a  unidirectional- 
fiber-reinforced  composite,  one  axis  is  obviously  the  fiber  direction. 


ar.G  Q  .•  j )  can  be  related  to 
_  , 

E^,  roisson  s  ratios 


If  one  wishes  to  express  the  stresses,  strains,  or  elastic  coefficients 
with  respect  to  an  orthogonal  coordinate  system  oriented  with  its  z  axis  in  line 
with  direction  X3,  with  its  x,y  axes  oriented  at  an  angle  3  with  respect  to  the 
X;,  :<2  axes,  standard  tensor  relations  can  be  used  to  obtain  appropriate  trans¬ 
formation  equations  i.1,2,3].  Then 


lV 


[V-;Y 


(7) 


'he  bar  denotes  that  the  quantity  is  related  to  the  (x,y,z)  coordinate  system. 
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Ic  is  noced  chat  in  the  array  of  equation  (7),  certain  of  the  coefficients 
(Cl6»  C26 »  and  0^5)  that  were  identically  zero  in  the  orthocropic  case,  equation 
(3),  are  no  longer  zero.  Thus, 


sll  S12 

§12  Szz 

§13  §23 

0  0 

0  0 

Sis  §25 


51  3  0 

52  3  0 

53  3  0 

0  S  4  4 
0  S45 

§36  0 


0  s16 

0  §26 

0  §36 

S45  0 

S55  0 

0  §60 


(8) 


The  number  of  non-zero  coefficients  has  increased  from  12  to  20  and  the  number 
of  independent  ones  has  increased  from  9  to  13.  The  elastic  symmetry  class  re¬ 
presented  by  equation  (8)  is  called  monoclinic.  In  ocher  words,  a  macerial  that 
is  orthotropic  with  respect  to  symmetry  axes  1,2,3  behaves  like  a  monoclinic 
macerial  with  respect  to  axes  x,y,z  as  defined  above. 


Table  2  summarizes  the  various  symmetry  classes  of  importance  in  fiber- 
reinforced  composite  materials. 


Table  2.  Summary  of  Selected  Elastic  Symmetry  Classes 


Symmetry  Total  No.  of  Non-Zero 

Class  Elastic  Coefficients 

Three-Dimensional-Stress  Cases: 


General  anisotropic  36 

Monoclinic  '  20 

Orthocropic  12 

Transversely  isotropic  12 

Isocropic  12 

Plane-Stress  Cases: 

General  anisotropic  = 

aonoclinic  (symm.  about  x 3=0)  9 

Orthotropic  (=  transversely 

isotropic  in  X2X3  plane)  5 

Transversely  isocropic  in 

xix2  plane  (=  isocropic)  5 


No.  of  Independent 
Elastic  Coefficients 


21 

13 

9 

5 

2 
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2.2  Stiffness  of  Laminates 

To  obtain  certain  desired  combinations  of  properties,  such  as  stiffnesses 
and  strengths,  ic  is  usually  necessary  to  arrange  individual  layers  (sometimes 
called  laminae)  of  fiber-reinforced  composite  materials  to  form  a  multi-layer 
laminate.  It  is  conceptually  possible  to  treat  structural  elements  constructed 
of  laminates  by  applying  the  relations  of  chree-dimensionai  elastodynamics  to 
each  individual  layer  in  conjunction  with  appropriate  continuity  and  boundary 
conditions.  However,  the  resulting  analysis  is  conceptually  so  complicated  that 
it  is  not  practical  to  carry  out  except  in  a  few  very  special  cases.  To  avoid 
this  complexity,  ic  is  the  usual  engineering  practice  to  make  certain  hypotheses 
regarding  the  kinematics  of  deformation  and  then  to  work  with  stiffnesses  which 
are  integrated  through  the  thickness. 
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The  simplest  kinematic  assumption  is  the  Bernoulli-Euler  one  for  slender 
beams  (or  the  Kirchhof f-Love  one  for  thin  plates  or  shells) :  plane  sections  re¬ 
main  plane  and  normal  to  the  deflected  middle  surface  during  deformation,  and 
undergo  no  thickness  stretching.  In  view  of  the  low  shear  modulus  of  most  com¬ 
posite  materials  relative  to  their  in-plane  elastic  moduli,  it  is  more  accurate 
to  use  the  Bresse-Timoshenko  hypothesis  [4,5]  for  beams  (or  the  Reissner-Mindlin 
one  [6,7]  for  plates  and  shells).  In  this  hypothesis,  it  is  assumed  chat  plane 
sections  remain  plane  but  not  necessarily  normal  to  the  deflected  middle  surface 
and  suffer  no  thickness  stretching.  Thus,  the  additional  deflection  due  to 
chickness-shear  deformation  is  accounted  for  in  an  approximate  sense.  Mathe¬ 
matically 


u(x,y,z)  =  uQ(x,y)  +  zB^x.y) 
v(x,y, z)  =  vQ(x,y)  +  z8v(x,y) 


w(x,y , z)  =  w  (x,y) 
o 


(9) 


Here  u,v  are  displacements  in  the  x,y  directions;  w  is  the  normal  deflection; 
uo,v0>wo  are  c^e  ®iddle-surf ace  displacements  in  the  x,y,z  directions;  (  ),  = 

?(  )/3x;  3X  and  3V  are  the  bending  slopes  and  differ  from  the  total  slopes  by 

the  respective  shear  strains.  Using  equations  (9)  in  the  equations  of  linear 
elasticity  theory,  one  gets 


-ti(x,y,z)  »  e“(x,y)  +  ZK.(x,y)  (i=l,2,6) 


(10) 


Here  s?  are  che  mid-plane  stretching  and  shearing  strains  and  < .  are  the  bending 
and  twisting  curvatures. 


Now  the  stress  resultants  and  stress  couples  are  related  to  the  correspond¬ 
ing  stress  components  as  follows  for  a  flat  laminate  (plate)  : 


/n/2 

N.  »  !  5.  dz  (i-1,2,6) 

1  * -h/2  1 


th/2 

Q,  =  ;  0.  dz  (i-4,5) 

J-h/2  1 


rh/2 

M  -  !  dz  (i-1,2,6) 

1  J -hi  2  1 


(11) 


Here  N,  are  in-plane  forces  per  unit  width,  Q^_  are  thickness  shear  forces  per 
unit  width,  M^  are  bending  or  twisting  moments  per  unit  width,  and  h  is  the 
total  laminate  thickness. 


Substituting  equations  (5)  and  (10)  in  equations  (11)  and  performing  the 
integrations,  one  obtains  the  following  laminated  plate  constitutive  relations: 


f  ^ 

A, 

'  3 

1  i 

•  = 

ij 

■  ij 

1  — 

1 

_Bij 

lDij 

1  — 

{Qi}  "  CA±  3(s  }  U,j-4,5) 


(12) 


Here  the  A^j  (i,j*l,2,6)  are  in-plane  (stretching  and  shearing)  stiffnesses,  A,-  4 
(i,j»4,5)  are  thickness-shear  stiffnesses,  3^  are  in-plane/out-of-plane  coupling 
stiffnesses,  and  Djj  are  out-of-plane  (bending  and  twisting)  stiffnesses  given 
by: 
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(13) 


(i, j=4.5) 


Here  index  l  refers  to  che  l-th  layer,  n  =  total  number  of  layers,  z._1  and  z. 
are  z-coordinace  positions  of  the  top  and  bottom  surfaces  of  the  typical  layer 
i,  and  the  kj_  are  the  shear-correction  coefficients  to  account  for  the  nonlinear 
distribution  of  the  thickness  shear  strains  through  the  total  thickness. 

The  basic  difference  between  the  theories  of  Refs.  [8j  and  [9]  is  that  the 
former  defined  the  plate  stiffnesses  (A^j ,  etc.)  on  the  basis  of  the  three- 
dimensional  material  stiffness  coefficients  (C^)  while  the  latter  defined  them 
in  terms  of  the  Qij  as  in  equations  (13).  Wang  and  Chou  [lO]  showed  that  the 
latter  approach  gives  more  accurate  results  even  for  relatively  thick  plates. 

There  is  a  nearly  unlimited  number  of  combinations  of  layer  orientations 
possible  in  a  laminate.  Table  3  lists  commonly  used  lamination  arrangements 
consisting  of  multiple  plies  all  having  the  same  thickness  and  of  identical 
composite  material  (except  for  orientation) . 


Table  3.  Characteristics  of  Some  Commonly  Used  Lamination  Arrangements 


Laminate  Name 

OescriDtion 

Aligned  parallel-ply 

All  layers  at  0° 

Off-axis  parallel-ply 

All  layers  ac  5 

Symmetric  cross-ply 

Alternating  0°  4  90° 
layers;  n=oc 

Antisymmetric 

cross-ply 

Alternating  0°  St  90° 
layers;  n=even 

Symmetric  regular 
angle-ply 

Alternating  +9  4  -6 
layers;  n=odd 

Antisymmetric  regu¬ 
lar  angle-ply 

Alternating  +9  4  -5 
layers;  n=even 

Symmetric  balanced 
angle-ply 

Canter  2  plies  at  5; 
remainder  alter¬ 
nating  +9  4-9; 
n*even 

Quasi- iso tropic 
(isotropic  A.,  only) 

ij 

One  or  more  sets  of 
layers  at  0,9;, 29;, 
3c; , . . , (X-l) 6; ,  where 
9;-lS0'J/M,  N  i  no.  of 
different  orienta¬ 
tions  (integer  _>  3) 

Stiffness  Characteristics 

Ref. 

A;e=A25=D;5=D'5=all  3^=0 

[11] 

All  31j=«0 

[11] 

Ai ;=A-;=D- ;=D-s*all  3..=0 

[3,11, 

J 

12] 

A,s»A-;-D^-D-c»0;  3-  -  —3  ; 

[3,11, 

remainder  of  3,- 

12] 

All  3^-0;  A;e  ,A;=  ,0;6  ,D;; 

[3,11, 

decrease  as  n  is  increased 

12] 

Al5*A25“B;;»B;;-3;;-B5g» 

[3,11, 

d1S“dL6“° 

12] 

A;;*A2s=all  B^j«0;  D;;  and 

[13] 

decrease  as  n  is 
increased 


A;:=A;--A,  Al:»vA, 

Agj ■(!-'-')  (A/2)  ,  A15-A;s-0 


[3,11] 
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2.3  Stiffness  of  Composites  with  Different  Behavior  in  Tension  and  Compression 


Analysis  of  materials  having  different  stress-strain  relations  in  tension 
and  compression  dates  back  to  Saint-Venant  (1864);  see  [l4].  Examples  of  mate¬ 
rials  which  have  been  found  to  demonstrate  this  behavior  include  various  cords 
and  fabric  embedded  in  rubber,  sintered  metal,  and  certain  biological  materials 

[14] ,  Modern  treatment  of  the  topic  daces  from  the  1965  work  of  Ambartsumyan 

[15] ,  in  which  he  instituted  the  concept  of  a  bimodulus  model  for  isotropic 
materials  with  different  properties  in  tension  and  compression.  This  concepc 
was  later  generalized  to  anisotropic  material  by  the  present  investigator  [16, 
17]  using  a  fiber-governed  symmetric  elastic  model  which  was  found  to  agree  well 
with  experimental  results  for  cord/rubber. 


In  Che  fiber-governed  symmetric  elastic  model  [16,17],  it  is  assumed  that 
there  is  a  complete  array  of  stiffnesses  which  govern  when  the  f iber-direction 
strain  (sf)  is  tensile  and  another  complete  array  when  it  is  compressive.  That 
is,  for  a  thin  laminate,  the  plane-stress  reduced  stiffness  matrix  is 


- 


ts.ji1  if  -  0 
feljzl  “  '£  <  ° 


(14) 


Here  the  index  k  refers  to  either  the  tensile  (k=l)  or  compressive  (k=2)  sec  of 
stiffnesses. 


For  a  beam,  place,  or  shell  undergoing  bending,  the  locus  of  points  at 
which  -tj=0  is  called  the  neutral  surface.  In  an  ordinary  (not  bimodulus)  mate¬ 
rial,  the  neutral-surface  location  is  not  especially  important.  However,  in  a 
bimodulus  macerial,  its  location  is  crucial,  and  in  general  it  does  not  coincide 
wich  the  middle  surface  for  such  materials.  Thus,  even  a  single  layer  of  bi¬ 
modulus  material  exhibits  bending-stretching  coupling  (analogous  to  an  unsym- 
mecrically  laminated  laminate) ,  except  when  the  magnitude  of  the  middle-surface 
strain  in  the  fiber  direction  exceeds  the  bending  strain.  In  the  latter  excep¬ 
tion,  a  bimodulus  material  acts  no  differently  than  an  ordinary  one. 


Depending  upon  the  structural  configuration,  boundary  conditions,  and 
loading,  even  in  a  single-layer  plate,  the  neutral  surface  in  general  is  not  a 
horizontal  plane.  Since  the  layer  properties  (Q^jk)  depend  upon  the  sign  of 
(and  thus  upon  the  neutral-surface  position)  and  the  plate  stiffnesses  (A,-i, 
3ii,Dij)  uPon  Qijk  as  an  equations  (13),  it  is  clear  that  in  general  a  bi- 
modulus  plate  is  nonnomogeneous  in  its  plane,  i.e.,  the  plate  stiffnesses  de¬ 
pend  upon  position  (x,y)  in  the  plane  [lS].  However,  unlike  a  plate  with  taper¬ 
ing  thickness,  here  the  plate  stiffnesses  are  not  even  known  a  priori.  This 
situation  is  reminiscent  of  the  elastoplas tic  bending  problem,  for  which  the 
elascoplastic  boundary  is  not  known  a  priori. 


2.4  Damping  of  Filamentary  Composites  and  Laminates 


Although  there  are  enumerable  mathematical  models  which  have  been  proposed 
to  describe  internal-friction  in  materials  [l9],  by  far  the  most  widely  used  is 
the  Kimball-Lovell  model  [20].  This  model  assumed  that  the  energy  dissipated  is 
proportional  to  the  square  of  the  amplitude  but  independent  of  frequency.  In  the 
case  of  sinusoidal  excitation,  it  is  convenient  to  apply  this  model  by  replacing 
the  elastic  stiffnesses  (Qj,j)  of  a  perfectly  elastic  material  by  the  complex 
stiffnesses  (Q^j)  defined  by 


Q 


ij 


(15) 
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Here  i  =  /-T,  QijR  =  real  (storage)  component  of  stiffness,  and  =  imaginary 

(loss)  component.  Although,  as  discussed  in  [2l],  there  are  considerable  dif¬ 
ficulties  in  applying  the  complex-stiffness  approach  to  nonsinusoidallv  excited 
systems,  no  alternative  proposed  to  date  appears  to  offer  much  promise,  vit'n  the 
possible  exception  of  the  fractional-derivative  approach  originated  by  Caputo  and 
recently  applied  to  isotropic  adhesive  material  by  Bagiev  and  Torvik  [22]. 
Preliminary  attempts  to  apply  this  model  to  boron-epoxy  have  not  been  success- 
ful. 


The  pioneering  work  in  analysis  of  damping  in  fiber-reinforced  composite 
materials  was  due  to  Hashin  [23],  using  the  elastic-viscoelastic  correspondence 
principle.  Unfortunately,  this  approach  requires  explicit  expressions  for  the 
elastic  stiffnesses,  and  such  expressions  are  generally  not  sufficiently  accu¬ 
rate  for  transverse  and  shear  stiffnesses  (Q22  and  Qss)-  These  limitations  were 
removed  in  [24],  a  complete  characterization  of  boron/epoxy. 

A  variety  of  experimental  techniques  have  been  used  to  experimentally 
characterize  the  dynamic  stiffness  and  damping  of  composite  materials.  Refer¬ 
ence  [25]  discussed  the  free  vibration,  pulse  propagation,  and  forced  vibratory 
response  methods.  Selected  recent  sources  of  experimental  data  for  various 
filamentary  composites  are  listed  in  Table  4. 

Table  4.  Selected  Sources  of  Experimental  Data  on  Damping 
of  Fiber-Reinforced  Composite  Materials 


Fiber  Material 

Fiber  Form 

Matrix  Material 

Ref. 

Glass 

Continuous 

Epoxy 

26,27 

3oron 

Continuous 

Epoxy 

23,29 

Graphite 

Continuous 

Epoxy 

28,30 

Glass 

Short 

Polyester 

31 

Aramid 

Woven  cloth 

Polyester 

32 

3.  BEAMS,  CURVED  3ARS ,  AMD  RIMGS 

Reference  [33]  surveyed  the  static  behavior  of  such  structural  elements 
and  discussed  25  pertinent  references  through  1973.  Here  emphasis  is  placed  on 
the  structural  dynamics  aspects. 

3. 1  Straight  Beams  of  Anisotropic  Material 

For  analyzing  the  lower  flexural  modes  of  straight  beams  made  of  ortho¬ 
tropic  material  (aligned  fiber-reinforced  composite  material),  the  dynamic 
version  of  the  simple  Bernoulli-Euler  type  theory  is  adequate.  Of  course,  the 
appropriate  Young's  modulus  (the  longitudinal  one)  must  be  used  in  defining  the 
flexural  rigidity. 

Since  composite  materials  typically  have  much  lower  ratios  of  shear  modulus 
to  Young's  modulus  than  homogeneous  isotropic  materials,  they  usually  exhibit 
more  transverse  shearing  action  than  isotropic  materials.  In  the  aligned  (.ortho¬ 
tropic)  case,  one  may  use  the  equations  of  3resse-Timoshenko  beam  theory  [4,5], 
provided  that  one  uses  the  the  appropriate  elastic  moduli:  longitudinal  Young's 
modulus  for  E  and  the  appropriate  (thickness  shear  or  in-plane  shear)  shear 
modulus  for  G.  Apparently  this  approach  was  used  for  orthotropic  beams  first  by 
Mowinski  [34]  and  later  by  others  l 35—36] . 

To  improve  the  dynamic  analysis  of  orthotropic  beams  still  further,  it  is, 
of  course,  necessary  to  use  a  three-dimensional  eiastodynamic  approach.  Such  a 


procedure  was  used,  for  example,  by  Ohnabe  and  Nowinski  [37],  who  extended  the 
Pochammer-Chree  equations  of  isotropic  elastodynamics  to  the  cylindr ically 
orthotropic  case. 

For  beams  made  of  off-axis  or  nor.alignea  fiber-reinforced  material  (i.e., 
general  anisotropic  material) ,  there  is  an  elastic  coupling  between  the  flex¬ 
ural  and  twisting  actions.  This  was  analyzed  by  Lekhnitskii  [2],  sections  27 
and  28,  for  the  static  case.  Apparently  the  first  dynamic  analysis  including 
this  coupling  was  Abarcar  and  Cunniff's  1972  analysis  [38],  who  combined  3resse 
Timoshenko  beam  theory  with  Lekhnitskii' s  flexural-torsional  coupling  theory. 
They  formulated  the  problem  as  one  in  Myklestad's  transfer-matrix  theory  [39]. 
Further  experimental  and  transfer-matrix  analyses  were  carried  out  by  Ritchie 
et  al.  [40-41],  Later  Miller  and  Adams  [42]  obtained  a  purely  analytical  solu¬ 
tion  for  a  theory  based  on  3ernoulli-Euler  beam  theory  rather  than  the  Bresse- 
Timoshenko  one. 

Recently  Teoh  and  Huang  [43]  presented  a  closed-form  solution  for  free 
vibrations  of  Bresse-Timoshenko-Lekhnitskii  beams  and  Teh  and  Huar.g  [44]  ana¬ 
lyzed  the  same  problem  using  finite  elements  based  on  an  energy  formulation. 
Most  recencly  Wallace  and  Bert  [45]  extended  the  transfer-matrix  approach  to 
include  Kinball-Lovell  material  damping.  At  a  considerable  saving  in  compu¬ 
tational  time,  they  were  able  to  predict  all  but  the  pure  "plate  modes"  and 
associated  frequencies  and  damping  factors  for  a  wide  beam  previously  investi¬ 
gated  experimentally  [46]  and  analytically  [47]. 

Mindlin  [48]  presented  a  three-dimensional  dynamic  analysis  of  the  low- 
frequency  modes  of  anisotropic  elastic  bars  of  arbitrary  cross  section  and 
applied  it  to  bars  with  elliptic,  triangular,  and  rectangular  cross  sections. 

3 . 2  Laminated  and  Sandwich  Straight  Beams 

A  laminated  beam  is  one  having  two  or  more  layers,  with  all  layers  of  the 
same  thickness  but  of  different  materials  (or  fiber  orientations).  In  contrast 
a  sandwich  beam  is  one  having  one  (or  more)  thick,  flexible  core  and  two  (or 
more)  thin,  stiff  facings.  A  simple  sandwich  has  one  core  between  two  facings, 
while  a  multicore  sandwich  ("club  sandwich")  has  more  than  one  core  sandwiched 
between  more  than  two  facings  [49].  Actually  a  multilayer  beam  with  constraine 
viscoelastic  layers  can  be  considered  to  be  a  multicore  sandwich. 

Although  laminated  beams  are  simpler  in  geometry  than  laminated  plates, 
relatively  few  vibrational  analyses  of  this  configuration  have  appeared  in  the 
recent  literature,  exemplified  by  [50-54], 

Since  there  is  a  voluminous  literature  on  the  linear  analysis  of  free  and 
sinusoidally  forced  vibration  of  simple  sandwich  beams  (much  of  the  early  work 
was  reviewed  in  [55])  ,  attention  here  is  focused  on  the  literature  of  the  past 
eight  years  with  particular  attention  to  various  complexities.  The  effects  of 
concentrated  masses  were  investigated  in  [56-57],  while  Ref.  [58]  considered  th 
effect  of  local  shearing  prevention,  such  as  that  occurring  due  to  the  presence 
of  rivets  near  the  ends.  The  effects  of  various  loadings  have  been  investigated 
these  included  impact  [59],  moving  force  systems  [60],  and  flutter  [61].  Hyer 
et  al.  [62]  investigated  the  effect  of  geometric  nonlinearity,  both  analytical! 
and  experimentally. 

The  effects  of  a  variety  of  configuration  parameters  on  sandwich-beam 
vibration  have  been  studied.  Sandman  [63]  considered  the  effects  of  segmenting 
i.e.,  discrete  changes  in  stiffness  and  mass  along  the  length  of  the  beam,  whil 
D.K.  Rao  and  Stuhler  [64]  investigated  the  effect  of  a  smoothly  tapering  cross 
section.  The  presence  of  pre-twist,  as  in  an  aerodynamic  surface  such  as  a 
compressor  blade,  was  shown  by  D.K.  Rao  [65]  to  result  in  drastic  reduction  of 
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che  effective  loss  faccor  for  Che  first  flexural  mode.  An  unsymmecric  sandwich 
is  one  which  does  noc  have  mid-plane  symmetry.  Vibrational  analyses  of  unsym- 
metric  sandwich  beams  having  different  top  and  bottom  facings  were  presented  by 
Yan  and  Dowell  [66],  Y.V.K.S.  Rao  and  Nakra  [67],  and  by  D.K.  Rao  [68-69]. 
Reference  [69]  also  included  the  effect  of  multiple  cores.  Y.V.K.S.  Rao  [ 7 0] 
analyzed  an  unsymmetric  sandwich  beam  that  he  called  a  dual-core  sandwich  con¬ 
figuration;  it  consisted  of  two  different  core  materials  (back-co-back)  sand¬ 
wiched  between  only  two  facings. 

3 . 3  Curved  Beams  and  Rings 

The  governing  equations  of  motion  for  curved  beams  and  rings  are  identi¬ 
cal;  che  difference  comes  in  via  che  boundary  conditions.  A  curved  beam  or 
arch  has  specified  boundary  conditions  at  each  end,  while  a  ring  must  have  a 
solution  which  is  circumferentially  continuous,  i.e.,  periodic  in  the  circum¬ 
ferential  angular  position.  Although  there  have  been  many  vibrational  analyses 
of  arches,  curved  beams,  and  rings,  apparently  che  only  one  strictly  applicable 
to  a  composite  ring  is  [7l].  It  treated  free  vibration  of  an  orthocropic- 
macarial  thick  ring  subjected  co  steady  centrifugal  loading,  supported  by  con¬ 
tinuous  radial  elastic  support,  and  included  transverse  shear  deformation.  The 
application  was  a  hoop-wound  rim-type  flywheel  for  automotive  energy-storage 
applications . 

Thera  have  been  a  few  dynamic  analyses  of  curved  sandwich  beams  and  rings. 
Ahmed  [72]  performed  a  finite-element  analysis  of  a  curved  beam.  References 
[73-75]  analyzed  the  steady-state  dynamic  response  of  damped  sandwich  rings, 
while  Sagartz  [76-77]  considered  the  transient  response. 


i.  FLAT  PANELS  (PLATES) 

Since  publication  of  Leissa's  comprehensive  monograph  [73]  on  free  vibra¬ 
tion  of  plates  in  general,  there  have  been  a  number  of  review  articles  [33,79-31j 
which  surveyed  the  vibration  of  composite-material  and  sandwich  plates.  Thus, 
in  che  interests  of  brevity,  only  a  few  very  recent  developments  not  discussed 
in  the  above  references  are  considered  here. 

■i .  1  Plates  of  Anisotropic  Material 

Linear  analyses  included  the  following:  for  orthotropic  rectangular 
plates,  Wilson  [82]  considered  che  response  to  moving  loads,  while  for  free 
vibration  of  anisotropic  rectangular  places  subjected  to  in-plane  forces,  Laura 
and  Luisoni  [83]  presented  a  Ritz-type  solution  with  a  polynomial  approximating 
function.  Reference  [84]  considered  moderately  thick  rectangular  plates  of  bi¬ 
modulus  orthotropic  material.  Reference  [85]  discussed  the  application  of  a 
combination  of  the  finite-strip  method  with  the  deflection-contour  method  and 
obtained  reasonably  good  results  for  an  orthotropic  square  plate. 
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latest  survey  [81j.  Rao  and  Singh  [92]  presented  an  optimal  design  synthesis 
procedure  to  design  rectangular  plates  of  minimum  weight  subject  to  contraints 
on  natural  frequencies.  A  finite-element  analysis  including  thickness-shear 
flexibility  was  developed  by  Reddy  [93]  and  compared  with  various  results  in 
the  literature.  Crawley  |_9d]  reported  extensive  experimental  results  for 
cantilever  rectangular  plates  of  graphite/epoxy  and  graphite/ epoxy- aluminum  and 
compared  the  results  with  moderately  thick  finite-element  predictions. 
Chatterjee  and  Kulkarni  [95]  analyzed  flutter  of  moderately  thick  rectangular 
panels  of  damped,  laminated  composite  materials. 

Elishakoff  and  Stavsky  [96]  analyzed  vibration  of  laminated  cylindrical!'/ 
orthotropic  circular  annular  places. 


5.  CYLINDRICAL!.!'  CURVED  PANELS  AND  COMPLETE  CYLINDERS 

This  area  has  been  surveyed  in  Leissa's  comprehensive  monograph  [97]  on 
free  vibration  of  shells  in  general  and  in  a  number  of  review  articles  [33,98, 
99]  on  vibration  of  composite-material  shells.  Hera,  in  the  interests  of 
brevity,  only  cylindrical  shells,  either  cylindrically  curved  panels  or  circum¬ 
ferentially  complece  cylinders,  are  reviewed.  Emphasis  is  placed  on  work 
carried  out  from  1973  to  the  present. 

It  should  be  recalled  that  for  places  there  are  primarily  two  classes  of 
analyses  (linear  and  nonlinear),  three  kinematic  situations  (3-D,  moderately 
thick,  or  chin),  and  two  material-geometry  cases  (single-layer/symmetrically 
laminated,  and  arbitrarily  laminated).  In  Che  case  of  shells,  there  is  a  much 
greater  variety  in  the  kinematic  situations.  The  hierarchy  of  shell  theories, 
from  most  accurate  to  least  accurate  (which  may  be  either  with  or  without 
thickness-shear  and  thickness-normal  action) [100]  is: 

1)  Exact  theory,  like  the  Langhaar-3oresi  theory 

2)  Second-approximation  theories,  like  Love's  seconc-approximacior. 
theory,  Flugge's  and  Novozhilov’s 

3)  First-approximation  theories,  such  as  Love's  first  approximation 
theory,  and  the  increasingly  popular  Sanders  theory 

4)  Morley's  shallow-shell  theory 

5)  Donnell's  very-shaliow-shell  theory 

The  selection  of  an  appropriate  theory  depends  upon  the  method  of  solution. 
For  example,  if  a  whole-shell  analytical  method  is  used  for  a  complete  shell, 
the  Donnell  and  Morley  theories  may  not  be  sufficiently  accurate  for  predicting 
the  frequencies  of  the  lower  circumferential  modes.  On  the  other  hand,  if  the 
shell  is  discretized  into  many  small  regions,  each  of  which  is  individually  very 
shallow,  then  the  Morley  or  Donnell  theory  may  be  entirely  adequate  and  the 
complexity  of  the  more  elaborate  theories  may  not  be  justified. 

5.1  Shells  of  Anisotropic Material 


Since  this  theory  is  relatively  simple,  especially  for  f irst-approximation 
or  simpler  theory  and  orthotropic  material,  recent  emphasis  has  been  on  compli¬ 
cating  effects. 

3abu  and  Reddy  [lOl]  presented  one  of  the  relatively  few  analyses  of 
cylindrical  panels.  Fortier  [102]  considered  the  effect  of  external  pressure  on 
cylindrical  panei  vibration.  Dym  [ 103 ]  investigated  the  effect  of  this  same 
loading  on  complete  cylinders,  while  Penzes  and  Kraus  [104]  also  added  the 
effects  of  torsion,  axial  force,  and  rotation  about  the  axis  (centrifugal  load¬ 
ing).  Kuptsov  [105]  considered  the  affect  of  centrifugal  loading  due  to  rota¬ 
tion  about  a  shell  diameter  at  one  end  of  the  shell.  The  effects  of  internal 
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irrocational  flow  [106-107],  external  flow  [l08],  and  internal  swirling  flow 
[109-110]  on  orthocropic  pipelines  have  been  analyzed,  as  have  the  effects  of 
random  boundary- layer  internal-pressure  fields  [ill]. 

In  one  of  the  relatively  few  analyses  of  nonlinear  vibration  of  a  shell  of 
nonlinear  anisotropic  material,  Shahinpoor  [112]  analyzed  a  thin  tube  of 
Ericksen-Rivlin  material.  The  large-amplitude  (geometrically  nonlinear)  vibra¬ 
tions  of  an  orthotropic  cylindrical  shell  with  axially  capering  thickness  was 
recently  investigated  by  Ramachandran  f 11 3] . 

Concerning  moderately  thick  orthocropic  shells,  the  work  of  War bur ton  and 
Soni  [ll4]  on  harmonically  excited  vibration,  including  material  damping,  is  one 
of  the  most  extensive.  Mangrum  and  Burns  [115]  considered  the  effect  of  a  dis¬ 
continuous  pressure  loading  moving  at  constant  axial  velocity.  Jain  [116] 
analyzed  the  vibration  of  a  vertical  cylinder  partially  filled  with  fluid,  while 
Shvecs  and  Marchuk  [117]  considered  a  pipe  with  an  internally  flowing  fluid. 

All  of  the  work  just  mencior.ed  was  limited  to  orthocropic  material.  One 
of  a  few  vibrational  investigations  of  anisotropic  (monoclinic)  cylindrical 
shells  of  moderately  thick  walls  was  recently  reported  in  [113].  The  analysis 
was  formulated  in  terms  of  a  higher-order  shell  theory  which  retained  some  of 
the  exact  kinematic  features  of  Che  Langhaar-3oresi  shell  theory  and  yet  in¬ 
cluded  thickness-shear  deformation  as  in  Reissner-Maghdi  theory.  The  predicted 
natural  frequencies  were  found  to  agree  reasonably  well  with  experimentally 
determined  ones  for  a  cylinder  of  unidirectional  material  oriented  at  30 
degrees  to  the  axis  of  the  cylinder. 

5 • 2  Laminated  and  Sandwich  Shells 

The  effects  of  several  different  boundary  conditions  on  free  vibrations  of 
laminated,  circular  cylindrical  shells  were  investigated  analytically  by  Abhat 
and  Wilcox  [119]  and  Fortier  and  Rossettos  [120 j.  Reference  [119]  is  particular¬ 
ly  interesting  in  that  it  introduced  a  new,  improved  method  of  reducing  the 
general  eighth-order  shell  frequency  determinant  to  a  fourth-order  one.  The  re¬ 
sults  obtained  using  this  method  are  considerably  more  accurate  than  those  ob¬ 
tained  by  Tu's  reduction  method  [l2l].  In  [119]  the  new  method  was  applied  to 
cross-ply  shells  with  either  clamped  or  freely  supported  ends.  In  [120],  four 
different  edge  conditions  were  considered  for  both  cross-ply  and  angle-ply  lami¬ 
nates.  In  [122],  four  different  boundary  conditions  were  investigated  experi¬ 
mentally. 

In  one  of  the  relatively  scarce  analyses  of  r.oncircuiar  laminated  cylindri¬ 
cal  shells,  Xoor  [123]  applied  the  so-called  multi-local  method,  which  is  a  vari¬ 
ation  of  the  Hermit ian  difference  technique.  Optimization  problems  relating  to 
vibration  of  laminated  anisotropic  circular  cylindrical  shells  were  investigated 
in  [124-125]. 

A  number  of  investigators  considered  the  effects  of  various  kinds  of  ex¬ 
ternal  loadings  on  free  vibration  of  laminated  cylindrical  shells.  References 
[l26]  and  [l27]  considered  external  pressure;  Paaovan  [l28]  included  pressure, 
axial  load,  torque,  and  centrifugal  and  Coriolis  forces. 

3erger  [129]  analyzed  the  vibrations  of  an  infinitely  long,  layered  ortho¬ 
tropic  cylindrical  shell  in  an  acoustic  medium.  Muggeridge  and  3uckley  [130] 
carried  out  analytical  and  experimental  investigations  on  symmetric  balanced 
angle-ply  cylinders  in  a  fluid. 

There  have  been  a  number  of  analyses  of  moderately  thick  circular  cylindri¬ 
cal  shells,  i.e.,  those  with  thickness-shear  flexibility  explicitly  included. 
Fortier  and  Rossettos  [l3l]  considered  shallow  cylindrical!"  curved  panels  of 


cross-ply  lamination.  Sinha  and  A.K.  Rath  [ 132 ]  considered  the  sane  geometry 
subject  to  free  vibration  and  buckling.  B.K.  Rath  and  Das  C 133 ]  analyzed  cross- 
ply  laminated  complete  cylindrical  shells  by  three  shell  theories:  a  refined 
moderately  thick-shell  theory,  "classical"  t'nin-shell  theory  (Love's  first 
approximation),  and  Donnell's  thin-shell  theory. 

Sun  and  Whitney  [134]  presented  a  theory  for  moderately  thick  laminated 
cylindrical  shells.  This  theory  was  applied  by  C.T.  Sun  and  P.W.  Sun  [135-136] 
to  several  cases  of  suddenly  applied  loading,  using  classical  separation  of 
variables  and  the  Mindlin-Goodman  technique. 

Various  applications  of  the  finite- element  technique  have  been  made  to 
moderately  thick,  laminated,  circular  cylindrical  shells.  Shivakumar  and  Murty 
[137]  developed  a  ring-type  element  with  sixteen  degrees  of  freedom.  Crawley 
[94]  compared  FEM  and  experimental  results.  Bradford  and  Dong  [138]  devised  a 
refined  element,  in  which  a  number  of  elements  through  the  thickness  comprise  a 
laminate, and  applied  it  to  shells  under  initial  stress. 

There  have  been  relatively  few  three-dimensional  elastodynamic  analyses  of 
laminated  cylinders.  Two  examples  of  such  analyses  are  those  of  Srinivas  [139] 
and  Mulholland  and  Gupta  [140]. 

Recently  the  activity  in  vibration  of  sandwich  shells  has  diminished  in 
comparison  to  its  vitality  in  the  1960's.  However,  one  should  mention  the 
theoretical  and  experimental  investigations  by  Harare  and  Sandman  [141-142]  on 
sandwich  shells  with  graphite/epoxy  facings  and  the  nonlinear  parametric  vibra¬ 
tion  analysis  of  Popov  et  al.  [143]. 


6.  SUGGESTIONS  FOR  FUTURE  RESEARCH 

The  author  believes  that  the  following  needs  are  most  pressing  in  the 
areas  covered  by  this  survey. 

1.  More  realistic  mathematical  modeling  of  material  behavior  is  urgently 
needed.  This  includes  such  effects  as  material  damping,  stress-strain  non¬ 
linearity  in  shear,  different  behavior  in  tension  and  compression,  and  effects 
of  temperature,  humidity,  and  material  damage  on  the  stiffness  properties. 
Experimental  verification  of  the  improved  models  is  also  of  importance. 

2.  Means  for  more  realistic  incorporation  of  practically  important  local¬ 
ized  discontinuities  should  be  developed.  These  factors  include  local  doublers 
and  edge  reinf orcements,  access  ports  and  other  cutouts,  and  attached  localized 
masses . 


3.  A  comprehensive  and  comparative  assessment  of  the  numerous  laminated 
shell  theories  presently  available  is  sorely  needed.  The  goal  should  be  to 
determine  the  simplest  theory  necessary  for  practical  engineering  calculation 
of  the  structural  dynamic  response  of  composite  shell  structures. 

4.  Development  of  reasonably  general  design  data  to  guide  the  structural 
designer  in  application  of  composites  in  dynamically  loaded  structures  is  most 
urgently  needed. 
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